Staggered Dynamics in Antiferromagnets by Collective Coordinates 
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Antiferromagnets can be used to store and manipulate spin information, but the coupled dynamics 
of the staggered field and the magnetization are very complex. We present a theory which is 
conceptually much simpler and which uses collective coordinates to describe staggered field dynamics 
in antiferromagnetic textures. The theory includes effects from dissipation, external magnetic fields, 
as well as reactive and dissipative current-induced torques. We conclude that, at low frequencies and 
amplitudes, currents induce collective motion by means of dissipative rather than reactive torques. 
The dynamics of a one-dimensional domain wall, pinned at 90° at its ends, are described as a driven 
harmonic oscillator with a natural frequency inversely proportional to the length of the texture. 

PACS numbers: 75.78.Fg, 75.50.Ee, 85.75.-d 



New developments have created opportunities for 
using antiferromagnets (AFMs) as active compo- 
nents in spintronics devices [1]. AFMs are ordered 
spin systems which lack a macroscopic magnetiza- 
tion in equilibrium because neighboring spins com- 
pensate each other. Despite this striking differ- 
ence from ferromagnets, AFMs have many proper- 
ties that are analagous to those of ferromagnets. 
The realization of high-temperature antiferromag- 
netic semiconductors [2], the feasibility of engineer- 
ing antiferromagnetic domain walls [3], the demon- 
stration of tunneling anisotropic magnetoresistance 
(AMR) [4], the possibility for ultra-fast spin switch- 
ing by light in AFMs [5], and the absence of stray 
fields make AFMs interesting for spintronics appli- 
cations. 

In magnetic materials, currents induce torques 
on the magnetic moments [6]. In ferromagnets, 
these torques can be used to switch the magneti- 
zation, induce steady state precession in magnetic 
oscillator circuits, or move domain walls. The- 
oretical [7] and experimental [8] results indicate 
that current-induced torque effects are present in 
AFMs as well, and that these effects are of the 
same order of magnitude as in ferromagnets. How- 
ever, several aspects are fundamentally different. 
For instance, the dynamics in antiferromagnets are 
described by coupled equations of the staggered 
field and the (out-of-equilibrium) magnetization. 
Current-induced torques affect these variables dif- 
ferently. 

In AFMs, the staggered field may spatially vary 
and is influenced by external magnetic fields and cur- 
rents. Traditionally, understanding the complex be- 
havior of the temporal- and spatial-dependent order 
parameter requires solving a set of coupled equa- 



tions with many degrees of freedom. In this Let- 
ter, we formulate a conceptually simpler theory of 
how external forces influence the staggered field and 
magnetization dynamics in AFMs in terms of a few 
collective coordinates. Our description is based on 
the phenomenological theory of insulating AFMs [9] , 
extended to account for charge current flow [10]. It 
includes the effects of dissipation, external magnetic 
fields, and both reactive (adiabatic) and dissipative 
(non-adiabatic) current-induced torques in slowly 
varying inhomogeneous antiferromagnetic textures. 

Consider a basic AFM lattice consisting of 
two magnetic sublattices, with magnetic moments 
mi(r,t) and m2(r, t), so that the total magnetiza- 
tion is m(r,t) = mi(r, t) + m2(r, t), and the anti- 
ferromagnetic order parameter is l(r, t) = mi(r, t) — 
m2(r,t). In the absence of magnetic fields and tex- 
tures, the equilibrium magnetization vanishes and 
l(r, t) is finite and homogeneous. Below, we con- 
sider the dynamics of the magnetization vector and 
the unit Neel vector n(r, t) = l(r, t)//(r, t). 

To the lowest order in textures and magnetiza- 
tions, the AFM free energy reads [9] 
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where a and A are the homogeneous and inhomo- 
geneous exchange constants, respectively. H repre- 
sents the external magnetic field. From the free en- 
ergy (1) and the constraints |n| = 1 and m • n = 0, 
which are valid for temperatures well below the Neel 
temperature, we can construct the effective fields 
fn = —SU/Sn = An x (V^n x n) — m(H • n) and 
fm = —SU/Sm = —am + n x (H x n). In all our 
results, we may generalize the free energy (1) by 
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adding anisotropy terms, e.g., easy-axis anisotropy 

Hals et. al. [10] introduced phenomenological re- 
active (adiabatic) and dissipative (non-adiabatic) 
current-induced torque terms, as well as dissipation. 
With these additional terms, the equations of mo- 
tion are 

h = (7f^ - Giih) X n + 7/7(J • V)n, (2) 
m = (7f,-G2n + /37(J-V)n) xn + T,,, (3) 

where 7 is the gyromagnetic ratio, Gi and G2 
are phenomenological Gilbert damping parameters, 
and T] {j3) parametrize the adiabatic (non-adiabatic) 
current-induced torque terms. Throughout this pa- 
per, we disregard all non-linear terms that are con- 
tained in Tni [10]. Eqs. (2) and (3) are the AFM 
analogs to the Landau-Lifshitz-Gilbert-Slonczweski 
equation for ferromagnets. By combining these 
equations, the magnetization can be expressed in 
terms of the AFM order parameter, giving a closed 
equation for the staggered field vector n to the linear 
order in the out-of-equilibrium deviations m, 9tn, J, 
and H: 

^ = -n X H + Gifn + (7? + • V)n 

7 

+a[7fn-G2n + 7/^(J- V)n]. (4) 

Here 7 = 7/(1 + G1G2) is a modified effective gyro- 
magnetic ratio in the presence of dissipation. Eq. (4) 
is the starting point for deriving the collective coor- 
dinate equation of motion for AFMs. 

In ferromagnets, magnetic textures are often rigid, 
so that only a few, soft modes dominate the magne- 
tization dynamics, as in the seminal work of Schryer 
and Walker on domain wall motion [11]. The evolu- 
tion of these soft modes can be described by a finite 
set of collective coordinates. This approach greatly 
simplifies the understanding of complex magnetiza- 
tion dynamics, making it possible to approximately 
describe the dynamics at low energies by considering 
only a few soft modes. 

The collective coordinate approach has recently 
been applied to magnetization dynamics in ferro- 
magnets [12]. We now present how the equations 
of motion for the collective coordinates can be con- 
structed for AFMs. We transform Eq. (4) by re- 
quiring the time dependence of the Neel field to be 
described by a set of collective coordinates {hiit)}: 
n(r, t) = n(r, The time derivative of the 

staggered field is then h = bidb^n. Similarly, h = 
bidi).n-\-0{b'^), where the second term is disregarded 
in our linear response analysis since it is quadratic 
in the driving forces. 



The dissipation is described in Eq. (4) via the 
terms Gif^ and aG2n. The first term scales as 
GiA/(A^r), where A and r are characteristic length 
and time scales of the staggered field texture. The 
second term scales as aG2/T. In analyzing the rel- 
ative strengths of these dissipative terms, we use 
the fact that the homogeneous and the inhomoge- 
neous exchange constants are related through a ~ 
A/{l'^(P) [13], where d is the lattice constant and we 
have introduced the AFM order parameter 1 above. 
Dissipation in metallic ferromagnets is small since it 
arises from the spin-orbit interaction in combination 
with electron scattering [14]. It is likely that similar 
mechanisms in antiferromagnets are also weak, and 
that they have comparable effects on the staggered 
field and the magnetization: Gil ~ G2// <C 1. From 
this we can conclude that 7 7 and that the second 
dissipative term, aG2n, dominates in realistic sys- 
tems, where the typical size of the texture A is such 
that X ^ d. Hence Gif^ can be safely disregarded 
in the equation of motion (4). 

Our main result is the equations of motion for the 
soft modes: 

M'^ (bj ^ jaG2bj) = F\ (5) 

This equation is derived by introducing the collec- 
tive coordinates to Eq. (4), taking the scalar product 
with dbjii.^ and integrating over the space. The dy- 
namics are equivalent to the classical motion of a 
massive particle subject to dissipation-induced fric- 
tion and external forces. This equation is model in- 
dependent and can be used to determine the param- 
eters of antiferromagnets, e.g. the Gilbert damping 
G2 and the homogeneous exhange constant a, which 
are usually difficult to identify in experiments. In 
Eq. (5), M'^ = J dV{db^n ■ db^n) is the effective 
mass arising from the exchange interaction between 
the spins. 

The total force inducing motion of the collective 
coordinates, = -\-Fj-\-Fjj^ is a sum of the ex- 
change force Fj^ = aj'^ J dV{dbin • f^), the current- 
induced force Fj = jJ (iV[7a/3 (95.11 • (J • V)n) + 
{v + Gif3){dbin • (J • V)n)], and the external field 
force Fjj = 7/ dVti • (n x db^n). More generally, 
Fj^ can also be expressed as = aj'^db^U^ to in- 
clude the effective material-specific forces which act 
on the antiferromagnet through the exchange inter- 
action and magnetic anisotropy. F} includes the re- 
active and dissipative current-induced forces, both 
of which are important for the dynamics of the col- 
lective coordinates bi. Fjj represents the response to 
an external magnetic field. Note that in the linear 
response regime, only time varying external mag- 
netic fields affect the dynamics of the collective co- 
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ordinates in AFMs, in contrast to the situation for 
ferromagnets [12], making the cohective motion in 
AFMs more resistant to stray fields. 

We now apply the general collective coordinate 
description, Eq. (5), to an isotropic one-dimensional 
antiferromagnetic texture, an orientational domain 
wall [15], in which the antiferromagnet is pinned in 
the X and z directions at z = and z = A, re- 
spectively. The staggered field n(z,t) varies slowly 
in the z direction, see Fig. 1. The pinning can 
be achieved by placing the antiferromagnet in con- 
tact with ferromagnets, as schematically shown in 
Fig. 1(c). In general, the staggered field can be ex- 
pressed in terms of two angles 6 and (j): n{z^t) = 
{cos cos ^, cos 6 sin ^, sin 6}. In equilibrium = O^q 
and (/) = 0, with Oeq{z) = 7tz/{2X). Without loss 
of generality, we assume that the out-of-plane angle 
(j) remains zero when a current passes through the 
system, which gives an AFM texture varying in the 
x-z plane only. 



a) 



b) 



c) 
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FIG. 1. (color online) A one-dimensional antiferromag- 
netic texture pinned at a relative angle of 90° in the left 
and right reservoirs, (a) shows the equilibrium orienta- 
tion of the staggered field, (b) depicts how a current J 
exerts a torque on the staggered field vector, forcing the 
center coordinate zq to be displaced by r, and (c) shows 
schematically a setup of an antiferromagnet between two 
pinning ferromagnets. 

In the steady state regime with a constant current 
along the z direction, J = Ji, the solution of Eq. (4) 
is Os{z) = 1(1 - e^^)/(l - e^^), where Q = (3J/A. 
As a collective coordinate representing the softest 
mode, we use the deviation of the texture center, r, 
from its equilibrium position zq^ which is the point 
where the x component of the staggered field vector 
equals the z component, 0{r) = 7r/4. In equilibrium, 
when there are neither applied currents nor external 



fields, the center coordinate is zq = A/2. Motivated 
by the steady state solution ^5, expanding for small 
Q in the low current regime to the linear order in the 
deviation r from equilibrium, we use Eq. (5) with the 
ansatz that the staggered field can be fully described 
by the sine and cosine of a function 0{z^ r): 
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Using this ansatz and the equation of motion (5), we 
find that the deviation from equilibrium, r, obeys 



(7) 



where ujq = 7(10Aa)^/^/A is the natural frequency 
of the system, F = {a/20Ay^'^XG2 describes dissi- 
pation, and g is proportional to the external forces. 
There are two contributions to the external forces: 
One from the current, gj = 57[a7/3J+(7^+/3Gi) J]/4, 
and the other from time-varying external fields, 
gn = —5XHy/{27T). For DC currents, the reac- 
tive (adiabatic) force parametrized by 77 plays no 
role, and only the dissipative (non-adiabatic) force 
parametrized by /3 is important for the texture dy- 
namics. When the driving forces are independent 
of time, Eq. (7) describes damped harmonic os- 
cillations about a new perturbed position r^ew = 
—PJX'^/{SA). This solution is valid as long as 
rnew < A/2. Hence, using /3* = f3Jd/A = -0.005, 
the approach works well for systems with lengths up 
to several hundred lattice constants. 

Numerical values for the natural frequency can 
be estimated for AFM metals. For example, in 
FeMn, the inhomogeneous exchange constant is A = 
0.94- 10~^^ J/m [16], the lattice constant is (ipeMn = 
3.6 A [17], and the magnetic moment per sublattice 
is 1.65 /iB, with /Jjb being the Bohr magneton, giv- 
ing a natural frequency of approximately 1 GHz for a 
FeMn texture with a length of 100 lattice constants. 

In Fig. 2, the solution of the time-dependent equa- 
tion of motion for r, Eq. (7), has been compared to 
numerical results of a micromagnetic simulation of 
the coupled Eqs. (2) and (3), with the boundary 
conditions described in Fig. 1. The equations were 
first written in dimensionless form by scaling the z 
axis with the lattice constant (i, and the time axis 
with i = (7a/) Other dimensionless quantities, as 
well as the numerical values used in the simulation 
presented in Fig. 2, are summarized in Table I. 

Fig. 2 shows that the complex spatio-temporal 
dynamics of the antiferromagnetic texture can be 
described by the motion of the single soft mode r. 
Good agreement between the approximate collective 
coordinate approach and the numerical simulations 



4 




FIG. 2. (color online) Transient response of the deviation 
r from the equilibrium position zq = A/2, after a con- 
stant current has been applied at time t = 0. The anti- 
ferromagnetic texture shows damped oscillations around 
a new perturbed position. The magnitude of the per- 
turbation depends on the system length and the current 
density. 



is also achieved when the current density is doubled 
(not shown). Fitting the simple equation of motion, 
Eq. (7), to experimental data, e.g., from AMR mea- 
surements, can provide good estimates of the phe- 
nomenological parameters in AFMs. The generality 
of Eq. (5) also makes the collective coordinate ap- 
proach a powerful tool for investigating the dynam- 
ics of more complex AFM textures with more than 
a single soft mode. 

The staggered dynamics represented by the cen- 
ter coordinate r can be measured via the AMR 
effect. The magnitude of AMR in bulk antifer- 
romagnets is not known, but since the tunnelling 
AMR [4] is significant, we believe its bulk value will 
be too. A plausible assumption is that the sim- 
plest possible phenomenological model of antiferro- 
magnetic AMR is similar to AMR in ferromagnets, 
but with the AMR depending on the orientation of 
the staggered field rather than the magnetization: 
p(n) = po + Pam(n • i)^, where po is the isotropic 
resistivity and pani is the anisotropic resistivity. In- 
tegrating the resistivity over the system, using the 
ansatz in Eq. (6), and expanding to linear order in r, 
gives the effect of the AFM texture on the resistance 
as R{t) = + Pani [A/2 - 8r(t)/7r2]. Therefore, it 
should be possible to observe the effects of both DC 
and AC currents. For DC currents, the total resis- 
tance R{t) will be enhanced or reduced depending on 
the current direction. For AC currents, by sweeping 
the frequency, one should observe enhanced devia- 
tions of the resistance when the frequency equals the 
natural frequency of the texture. This setup offers 
the possibility of measuring the effect of the current- 



induced torque on the staggered field, a phenomenon 
which is, in general, difficult to observe experimen- 
tally. 

TABLE I. Dimensionless numerical constants 



Constant Composition Value 



* 

a 


al'^dP/A 


1 


ai 


Gil 


0.01 


0.2 


G2/I 


0.01 




pjd/A 


-0.005 


* 


r]Jd/{Al) 


-0.005 


A* 


X/d 


100 



We can also apply our collective coordinate ap- 
proach to an antiferromagnetic domain wall de- 
scribed by the Walker ansatz: tdiii{Ow) = e*^^"^^)/^^. 
Here A^^ = A/K^ is the domain wall width and 
Kz is the easy axis anisotropy. We also reintroduce 
the out-of-plane tilt angle and use the center of 
the domain wall r^^, ^^y, and the domain wall width 
as the three collective coordinates. In agreement 
with the simplified treatment in Refs. 10 and 18, by 
applying a constant current, the domain wall motion 
gradually relaxes to a steady state, where the wall 
moves with the constant velocity ~ — 7/3J/G2. 
Our approach shows that the out-of-plane tilt an- 
gle is a hard mode, which can only be excited by a 
time varying external magnetic field. This is very 
different from the motion of domain walls in ferro- 
magnets, where a moving domain wall also has a fi- 
nite tilt angle [19]. Additionally, in linear response, 
there is no distortion of the domain wall width for 
ant iferromagnet s . 

In conclusion, we have derived equations of mo- 
tions for the collective coordinates corresponding 
to soft modes of AFM textures to the linear order 
in currents, magnetization, and external magnetic 
field. In contrast to ferromagnets, the dynamics are 
second order in time derivatives, e.g., the effective 
particles described by the soft coordinates acquire 
a mass, and have no first-order contribution from 
time-independent external magnetic fields. We have 
applied our theory to a one-dimensional model of a 
slowly varying antiferromagnetic texture pinned at 
90° at the edges, and found the natural frequency 
and deviations of the center coordinate in terms of 
the system parameters. The results show that the 
dissipative (non-adiabatic) current-induced torque is 
crucial for the dynamics of the AFM textures. 
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